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, Abstract. In this note, first the Recurrent Variational Approach (RVA) is intro- 

' duced by using as example a non-trivial spin-model, the spin-1/2 antiferromagnetic 

two-leg-ladder. Then, a first application of this scheme to the electronic structure 
of conjugated polymers is proposed. An analogy between the two-leg ladder and 
the dimerized chain described by the Extended Peierls-Hubbard Hamiltonian is un- 
\ derlined and comparisons are made with DMRG calculations. The proposed ansatz 

^ . is useful to get some analytical insight and it is a good starting point for further 

CNJ ' improvements using the RVA. 
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1. Introduction 



Conjugated molecules have attracted the attention of chemists and 
\ physicists since the early days of quantum mechanics. A remarkable 

example is given by Hiickel who developed in the thirties the first 
Q , independent-electron theory applicable to polyenes (Hiickel, 1931). In 

the sixties, J. A. Pople and S.H. Walmsley (Pople and Walmsley, 1962), 
introduced by using this simple theory in their pioneering work what 
r> \ would become twenty years later the solitonic excitations. Following 

■ this line, Su, Schrieffer, Heeger and, independently, Rice, proposed a 

model specifically devoted to the polyacetylene - known as the SSH 
model - where the role of the electron-phonon interaction is especially 
emphasized (Su, 1979; Rice, 1979); here, the polyacetylene is thought 
of as a Peierls insulator and the solitons, polarons and bipolarons as 
the major relevant excitations in this system. A very nice review of 
these early studies is included in the book of L. Salem (Salem, 1966) 
and a review about the SSH model can be found in the next reference 
(Heeger, 1988). 

As a counterpart to this one-electron theory, the major role of 
electron correlation was soon recognized. In the fifties, the celebrated 
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Pariser-Parr-Popic (PPP) Hamiltonian was suggested as a possible 
candidate for the studies of conjugated oligomers (Pariser and Parr, 
1953; Pople, 1953). Later, Ovchinnikov et al. suggested to consider 
conjugated polymers as Mott-Hubbard instead of Peierls insulators 
(Ovchinnikov, 1973). Then, one has to face a N-body problem with 
a strong Coulomb interaction with a non-zero range; it was subject 
to many numerical works such as exact diagonalizations of very small 
oligomers (Chandross, 1999), selected CI calculations (Tavan and Schul- 
ten, 1987) and, quite recently, Density Matrix Renormalization Group 
studies (Bendazzoli, 1999; Boman and Bursill, 1998; Jeckelmann, 1998; 
Shuai, 1998). The PPP Hamiltonian was also a preferential model to 
support the development of new methods before applying them to more 
realistic models; a very good example of such work is given by the 
Coupled Cluster method first implemented in quantum chemistry for 
the PPP Hamiltonian (Cizek and Paldus, 1971). 

Today, it is well recognized that both the electron-phonon and 
electron-electron interactions are equally important (Baeriswyl, 1992). 
Moreover, interchain interaction and disorder are also of importance. 
Therefore, the problems one has to face to fully understand conjugated 
polymers stay quite challenging as it is the case for most of the recent 
organic compounds including high Tc superconductors. Despite a large 
amount of work, from the early thirties to now, the main questions 
remain awaiting an answer! For instance, a comprehensive picture for 
the ground state of these systems is still needed; it is the goal of this 
work to propose a relatively simple and, we believe, promising, way to 
understand the ground state of conjugated polymers. 

More generally, one dimensional correlated electron systems appear 
over the years as a growing and important field in condensed matter 
theory. One main reason of that comes from the fact that these systems 
serve as a theoretical laboratory to explore new methods of solution, 
analytically or numerically. A second main reason - we would say 
more physical - comes from the fact that over the years, more and 
more experimental realisations of these systems appear. New concepts 
emerge from such studies as the ones included in the phenomenology 
of Luttinger-liquids (Haldane, 1981). 

Among the technical methods proper to the one dimensional geome- 
try, one may cite the Bethe ansatz (Bethe, 1931), the bosonization tech- 
niques (Haldane, 1981), and, more recently, the Density-Matrix Renor- 
mahzation Group (DMRG) method (White, 1993; Peschel, 1999) and 
a closely related scheme which is directly considered in this note, the 
Recurrent Variational Approach (RVA) (Sierra and Martm-Delgado, 
1997; Peschel, 1999). The two first methods are analytical and the 
third one is numerical; the RVA method is in between. 
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The RVA, presented here, is a variational method invented recently 
in the context of spin-ladders (Sierra and Martm-Delgado, 1997). It 
is a method very closely related to the DMRG scheme with the main 
differences coming from the fact that only one state is retained as the 
best candidate for the ground state in the RVA on the contrary to the 
DMRG which considered much more states (Peschel, 1999). Most often, 
the results obtained with it are less accurate than the DMRG ones, but 
it is much easier to get a physical insight into the problem; the analysis 
of the results is simpler and the RVA could provide a "physical" picture 
of some phenomena which, wc believe, are of utmost importance in the 
field of conjugated polymers where the most important ingredient are 
still not fully recognized (Sariciftci, 1997). 

This note is organized as follows. In section II, the RVA is presented 
using the example of the two- leg spin ladders. In section III, the RVA 
method is applied to the study of the dimerized chain described by 
the Extended Peierls-Hubbard model. A nice similarity between the 
two-leg ladder and this system is pointed out and comparisons with 
DMRG calculations are made. We conclude in section IV and give some 
reasonable perspectives. 



2. RVA method applied to two— leg Spin Ladders 

A spin ladder is an array of coupled spin chains. The horizontal chains 
are called the legs, the vertical ones, rungs. In the case of spin one- 
half antiferromagnet spin-ladders, these systems show a remarkable 
behaviour in function of the number of leg: there is a gap in the excita- 
tion spectrum of even-leg ladders and, on the contrary, no gap in the 
excitation spectrum of odd-leg ladders. In terms of correlation lengths, 
this means that there is short (long) -range spin correlation in even 
(odd) -leg ladder (sec (Dagotto and Rice, 1996) for a review). 

The even-leg ladders show a spin-liquid ground state described by 
a Resonance Valence Bond (RVB) state introduced by P.W. Anderson 
(Anderson, 1987); such systems may be efficiently described by using 
a direct-space method as the DMRG and RVA ones. In this note, we 
will consider only the simplest even-leg ladder, the two-leg ladder, and 
we will describe it with the simplest RVB state, the Dimer-RVB state 
where the elementary singlets are between two nearest-neighbour sites. 
This problem is easily solved by using the RVA method (Sierra and 
Martfn-Delgado, 1997). Moreover, with the same scheme, it is possible 
straightforwardly to go beyond the simple Dimer-RVB state by consid- 
ering more extended elementary singlets (Sierra and Martfn-Delgado, 
1997), and to study more complex systems as the four-leg ladders for 
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instance (Roncaglia, 1999). In principle, this method should be rele- 
vant for every onc-dimcnsional system with exponentially decreasing 
correlation lengths as conjugated polymers (Pleutin, 1999). 

In this note we consider the simplest case given by the two-leg ladder 
(see figure (1)). The corresponding Hamiltonian is the following 

N-l N 

H = jJ2iS,. + • S^) + J'^ Sfe • (1) 

k=l k=l 

where J (J') are the exchange integral for the leg (rung) with J > 
(J' > 0). denotes the number of rungs of the ladder and open periodic 
boundary conditions have been assumed. 



1 J 2 


J c 


N- 


1 J N 


J' 


1 1 

J' 


1 1 

J' 


J' 



1 J 2 J 3 N-1 J N 

Figure 1. The two-leg ladder with an AF-Heisenberg model described by two 
coupling constants: J for the links in the legs, and J' for the vertical rungs 

In the limit where J' is much larger than J, the ground state of (1) 
is simply given by the product of singlets localized along the rungs (see 
fig. (3)). Starting from this limit, one may include fluctuations around 
it to treat the general case described by the Hamiltonian (1). A minimal 
way to do it, is to consider instead of singlet on rungs, pairs of nearest 
neighbour singlets on the leg; these local fluctuations are called resonon 
in the rest of the paper. Therefore, the resonance mechanism - in the 
sense of Pauling (Pauling, 1960) - we will consider in the following are 
represented in figure (2). 



< ► 



Figure 2. The basic bond resonance mechanism between horizontal and vertical 
bonds in an elementary plaquette of four sites. 

With the two selected local configurations, the singlet localized on 
a rung and the resonon, we have to generate the corresponding Hilbert 
space; the configurations are characterized by the number of resonon, 
M, and their positions, {xi}, as it is shown in the next examples. The 
Hilbert space is splitted according to the number of resonon; we have 
then to consider 
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■»>=1 I I I I I 

Figure 3. The zero-resonon state |0) in a two-leg ladder. It is made of all vertical 
rung singlets. It is the reference state in the strong coupling limit J' 2> J. 

— the zero-resonon sector with one unique state |0) (Fig. 3), 

— the one-resonon sector with the states (Fig. 4), 

'-=1 I I Z I I 

Figure 4- A generic one-resonon state in a two-leg ladder. 

— the two-resonon sector with the states {|2;i,X2)} (Fig. 5), 

= ! ~ I I ~ 1 I 

Figure 5. A generic two-resonon state |a;i,a::2) in a two-leg ladder, 
etc... 

— the Y~resonon sector with one unique state \xi,X2, ■ ■ ■ ,xm) (Fig. 6). 

With this classification of states we may write the following gener- 
alized Dimer-RVB state for a two-leg spin ladder, 

[N/2] 

\'>P{u))bvb = u'^^ ^ \xi,X2,...,Xm) (2) 

A/=0 {xi,X2,...,xm} 

where the sum is taken over all sectors with M resonons and u is 
the amplitude for a resonon. Here u is a variational parameter to be 
determined upon minimization of the ground-state energy. As a result, 
it is a function of the ratio of couplings, namely, u = u{J / J'). 

Since we kept only two very localized building blocks, it is possible 
to generate the state (2) in a recursive manner which we present now. 
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Figure 6. The all-resonon state |a;i, 2:2, . . . , a; jv ) in a two-leg ladder. 

2.1. Recurrence relations for the wave function 

Denoting by |A^) = \il){u))iivB a generic RVB state (2) for a two-leg 
ladder of N rungs, there are only two possibilities or movements to 
create generic RVB states of higher length, namely, 

— addition of one vertical rung to create the state + 1), 

— addition of one pair of horizontal bonds (resonon) to create |A^+2). 

From these arguments we can establish that the Dimer-RVB states 
(2) satisfy a recursion relation given by, 

\N + 2) = \N + 1)® \<j)l)N+2 +u\N)® \^2)n+i,n+2 (3) 

where the state denoted by |(/)i)iv+2 is a vertical rung at position N + 
2 and \(l)2) n+i,n+2 is made up of a pair of horizontal bonds located 
between the rungs at (A^ + 1, + 2), i.e.. 



|<^2)7V+i,iV+2 = (iV + 1, iV + 2)(iV + 2, iV + 1) (4) 

where [N, M) denotes a singlet between the sites N and M. 

Using (3) one can generate recursively the Dimer RVB state \N) 
from previous states with lower length and estimate variationally the 
ground-state energy as, 

- (2iV)(Ar|iV) -^^^^JIJ^^) (5) 
The value of u is fixed by minimization of (5). 

2.2. Recurrence relation for the norm 

To compute the norms of the states (2) let us define 



{N\N) (6) 



The two selected local configurations are not orthogonal, it is then 
convenient to define an auxiliary function as follows, 
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Yr, = {^(t>,\®{N-l\)\N) 

The second order RR's for the states leads (3) to a closed 
for the overlaps Z^, , namely, 

Zn+2 = Zn+i + uYn+1 + u^Zn 
where we have made use of the result, 

N+l{(l>l\4>2) N,N+1 = ^|0l)iV 

The RR's (8) together with the initial conditions, 

Zo = Zi = l , = 0, Fi = 1 (10) 

determine Z^ and Ijv for arbitrary values of as functions of u. 

Note that instead of introducing the function Yn, it is also possible 
to make first linear combination of our building blocks, making used of 
(9) (Eder, 1999), in order to guaranty the orthogonality between the 
selected local configurations. The overlap would then be given by only 
one recurrent relation as it will be the case in the next section. 

2.3. Recurrence relation for the expectation value of the 
Hamiltonian 

The expectation values of the Hamiltonian can be determined following 
the same steps outlined above. To this end, we introduce the quantities, 

En = {N\Hn\N) ,DN = U(t)i\®{N-l\)HN\N) (11) 

To obtain the RR's satisfied by E]\f,D]\i, one splits the Hamiltonian 
Hn of a two-leg ladder of length N (1) into two pieces: 

Hn = Hn-1 + Hn-i,n (12) 

where H^^i is the Hamiltonian of length N — 1 and Hn^i^n is the rest 
of the whole Hj^f Hamiltonian, i?Ar_i = Hi\f — ffjv_i, which is made 
of one vertical rung and two horizontal links. With this splitting and 
using (3) and (6), we obtain, 

En+2 = En+1 + J'eoZN+1 + u{Dn+i + (2 J + J')€oYn+i) 

W{EN + 2JeoZN) (13) 
Dn+2 = -E'at+i + J'eoZN+1 + + (2J + J')eoYN+i) 



(7) 

set of RR's 
(8) 

(9) 
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Table I. The values —e^/J' are mean field (MF) values taken 
from (Gopalan, 1994), while -e^^°/J' are Lanczos values taken 
from (Barnes, 1993). 



J/J' 


u 














0.375 


0.375000 




0.2 


0.128521 


0.383114 


0.382548 




0.4 


0.323211 


0.40835 


0.405430 




0.6 


0.578928 


0.44853 


0.442424 




0.8 


0.87441 


0.499295 


0.489552 




1 


1.18798 


0.556958 


0.542848 


0.578 


1.25 


1.58519 


0.63518 


0.614473 


0.6687 


1.66 


2.21853 


0.772172 


0.738360 


0.8333 


2.5 


3.39153 


1.06915 


1.002856 


1.18 


5 


5.9777 


1.99285 




2.265 



where eo = —3/4 is the lowest eigenvalue of the operator Si • 82- 
The initial conditions for Ej^ and Dj^ are 

Eo = 0, Ei = J'eo Do = 0, Di = J'eo- (14) 

The RR's for the energies involve the norms of the states and depend 
both on u and the coupling constants J', J. 

By working with orthogonal building blocks (Eder, 1999), it is also 
possible to go back to only one recurrent equation for the expectation 
value as for the norm; this will be the case in section 3. 

2.4. Results for the Variational Ground-State Energy 

The ground-state energy for a ladder of length N is estimated. In the 
thermodynamic limit N ^ 00 one can find a closed expression for the 
density energy per site (Sierra and Martm-Dclgado, 1997), 

- r 1 En _ R{a) 

~ nZo 2N Zn ~ 2aQ'{a)P{a) ^ ' 

in terms of three polynomials P,Q,R evaluated at the biggest root a of 
the cubic polynomial Q{y)- Finally one looks for the absolute minimum 
of Coo by varying the parameter u. In Table I we show the ground-state 
energies per site for different values of the coupling constant ratio J/ J' , 
varying through strong, intermediate and weak coupling regimes. 

In the strong coupling regime J/ J' < 1 the RVA states give a slightly 
better ground-state energy than the mean field result. This latter state 
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Table II. Here m is total the number of multiplets, Na is the number of independent 
variational parameters, dj is the number of multiplets with spin J, gMP.DMRG jg ^j^g g ^ 
energy density of the matrix-product state (DMRG), 1 — Pm is the probability of the states 

truncated out in the DMRG and (f^^ is the spin correlation length of the MP state. The 
exact results are given by Coo = 1.4014845 and ^ = 6.03 (White and Huse, 1993). 



m 


Na 


di/2 


d3/2 


•^6/2 


_ MP 


DMRG 
Coo 


1-Pm 


^MP 


1 





1 








1.333333 


1.333333 


1.6 X 10"^ 


0.910 


2 


2 


1 


1 





1.399659 


1.369077 


1.4 X 10"^ 


2.600 


3 


4 


2 


1 





1.401093 


1.392515 


1.3 X 10"^ 


3.338 


4 


7 


2 


2 





1.401380 


1.401380 


1.6 X 10"^ 


3.937 


5 


10 


2 


2 


1 


1.401443 


1.401436 


7.6 X 10"'^ 


4.085 


6 


13 


2 


3 


1 


1.401474 


1.401468 


1.3 X 10"" 


4.453 



produces rather unphysical results for J/ J' > 1, which does not occur 
in our case. The RVA results compare also well with the exact results 
for J' ~ J. 

The RR's can also be used to compute the spin correlator (Sj • Sj) 
which has an exponential decay behaviour exp(— |i — with ^ the 

spin correlation length which satisfies the equation 

^3^3 _ (2 + u)u^C'^ - (2 + 4u)u'^C + 4u^ = (16) 

where £. = e^/^. For J = J' one finds S_ = 0.737, which can be compared 
with its exact value given by 3.2. The latter results imply that at the 
isotropic point the bonds extended over three rungs are quite impor- 
tant. The simplest improvement of the Dimer-RVB state is to add a 
bond of length analogue to a knight move, shown in figure (7); this 
first improvement leads to a third order RR. The ground state energy 
per site of this state is given by -0.5713 while the spin correlation length 
is 0.959. Hence there is an improvement in both quantities but still one 
needs longer bonds. 




Figure 7. Local configurations extended over three rungs (the knight move) 

The RVA method has already been applied to several model sys- 
tems: to the Heisenberg two-leg ladder as we have seen above (Sierra 
and Martfn-Delgado, 1997) and to two generalizations of it, the t-J 
(Sierra, 1998) and the Hubbard (Kim, 1999) two-leg ladders where 
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the t-J and the Hubbard Hamiltonian, respectively, are defined on the 
lattice represented in the figure (1); it has been also applied to other 
ladders, the four-leg ladder (Roncaglia, 1999) and the diagonal ladders 
(Sierra, 1999), and to the spin-1 chain (Dukelsky, 1998). In this latter 
case, very accurate calculations have been made by considering sev- 
eral local configurations (instead of our two retained here); the results 
are compared with Density Matrix Renormalization Group (DMRG) 
calculations (White and Huse, 1993). We see that, in any case, the 
comparison are very good and even sometimes the RVA method gives 
better results than the DMRG (sec Table II). A more complete review 
about RVA can be found in the following reference (Peschel, 1999). 



3. RVA method applied to Conjugated Polymers 

The low-energy properties of conjugated polymers involve 7r-electrons. 
Usually, they are described by the use of effective models as the well- 
known Pariser-Parr-Pople Hamiltonian (Pariser and Parr, 1953; Pople, 
1953). This model includes both electron-phonon interaction (semi- 
classically) and a long range Coulomb interaction. The characteristic 
values for each kind of energetic contribution, kinetic, electron-phonon 
and electron-electron terms are approximatively of the same order of 
magnitude (Baeriswyl, 1992); in this regime, the study of the PPP 
Hamiltonian - or some other short versions of it - arc really not obvious. 

Compared to the previous model used in section (2), the spin 1/2 
antiferromagnetic Heisenberg model for a two-leg ladder, the PPP - or 
related - model is much more complicated: the long-range part of the 
Coulomb interaction is very difficult to treat and the dynamical variable 
are still not only spin but electron which involves spin and charge degree 
of freedom. However, due to the alternance between monomers, which 
can be large, and the link (the single bond) between the monomers, we 
will underline in this work a remarkable similarity between conjugated 
polymers and two-leg ladder. A similar wave function than (2) are 
proposed here for conjugated polymers - the rung of the ladder is played 
by the monomer and the resonon excitation of the two-leg ladder by 
a set of intermonomer fiuctuations (Pleutin and Fave, 1998; Pleutin, 
1999). We believe this wave function a natural and good starting point 
for further refinements. 

We choose as model system the simplest conjugated polymers - the 
trans-polyacetylene. This compound shows a dimerized structure with 
an alternance between double bond (1.35A) and single bond (1.45A); 
the monomer is then simply a double bond (see figure (8)). The vr elec- 
trons are assumed to be effectively describe not by the full PPP Hamil- 



paper.tex; 1/02/2008; 10:32; p. 10 



The RVA applied to the Electronic Structure of Conjugated Polymers 



11 



tonian but by the Extended Peierls-Hubbard model (EPH) (Baeriswyl, 
1992; Jeckelmann, 1998), a short version of it, for simphcity 

H = -^{t-{-'^T^){ci+l,a(^n,<j+(i,aCn+l,a) + U ^nn^nnl+V^nnUn+l 
n,a n n 

(17) 

The operator cli,a (c„ ^r) creates (annihilates) an electron of spin a at 

site n, Una- = Cn^aC^^ and n„ = n^t + n„j^. t is the nearest-neighbour 
hopping term without dimerization, A is the usual dimerization order 
parameter, U is the on-site Hubbard repulsion and V is the nearest- 
neighbour charge density - charge density interaction. We will note in 
the following for convenience, td = t+^ and tg = t — the hopping 
term for the double and the single bond respectively. 




Figure 8. Model for the trans-polyacetyleno and the system of bonding - antibond- 
ing levels associated with each double-bond 

Because of the dimerization, the one electron density is more pro- 
nounced on the double bonds than on the single bonds; this is more and 
more apparent when the dimerization increases. More generally, this is 
the case for every conjugated polymers where the one-electron density 
is peaked on the monomer region. In this sense, conjugated polymers 
are not strictly one dimensional systems but rather intermediate be- 
tween quasi-zero and quasi-one dimensional systems (Pleutin and Fave, 
1998; Soos, 1997; Mukhopadhyay, 1995). It seems then convenient to 
perform a unitary transformation which favours the orbitals localized 
on the monomers, in our case, on the double bonds. The new operators 
are then given by 



R(t) -J_fc(T) +c(t) \ /|(t)-J_fJt) -c^^) \ (m 

These operators create or annihilate vr electrons in the corresponding 
bonding or antibonding orbitals of the double bond n (see figure (8)). 

As for the case of the two-leg ladder, we select the Local Con- 
figurations (LC) which are the much important for the ground state 
wave function. This selection of the relevant LC is made by combining 



paper.tex; 1/02/2008; 10:32; p. 11 



12 



S. Pleutin, E. Jeckelmann, M.A Martm-Delgado, G. Sierra 



energetic and symmetry considerations. For that purpose, we use the 
electron-hole symmetry operator, J, to classify the LC. On the Fock 
space of a single site n, the action of this operator are summarized as 
follows (Ramasesha, 1996) 

4|0>=|n>, in I a >= (-in a >, 4|U>=-|0> (19) 

The electron-hole symmetry operators for several site is just given by 
the direct product of such single site operator J = j„. 

The electronic configurations can be splitted into two classes of 
electron-hole symmetry denoted by (+) and (-). The ground state is 
in the (+) sector, J \ GS >=| GS >; therefore we select the lowest LC 
in energy which belong to this class of symmetry. We kept only four 
of such LC, represented in figure (9) which were shown as the most 
important ones in finite cluster diagonalizations (Chandross, 1999). 
We present briefly now these LC; their associated creation operator 
and their energy. They are classified in two categories following their 
extension: the LC localized on one double bond are called Molecu- 
lar Configurations (MC) - to underline the analogy with molecular 
crystals; the LC extended over two nearest-neighbour double-bonds 
are named Nearest-Ncighbour-Intermonomer-Fluctuations (NNIF). In 
such classification, the singlet on the rungs of the previous section enters 
in the first class and the resonon in the second one. 

T T 

- ^ h + + 

^ ■ — • + ■ + + 

a b c d 

Figure 9. The set of local configurations considered in this work for the dimerized 
chain; a) the F-LC b) the _D-LC c) the Cti-LC and d) the TT-LC made by combining 
two localized triplet symbolized by the T on the figure. 



3.1. Molecular Configurations (MC) 

The MC are the LC completely localized on the monomers. They give 
an appropriated description of the system in the strong dimerized limit. 
There is two kinds of MC in the (+) electron-hole symmetry sector 
(Pleutin and Fave, 1998; Pleutin, 1999). 

— In the first, the monomer n is in its ground state; it is described 
by the following creation operator 
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4 = ^Lt^a (20) 

For a reasonable choice of parameters, this is the lowest LC in 
energy so that we choose as reference state 

I >= JJ Ft I Vacuum > (21) 



where | Vacuum > is the state without any tt electron. With 
respect to this reference state, the local operator Fn is simply given 
by the Unity /„. We have named this LC, F-LC. In the following 
all the creation operators and the energies are defined with respect 
to I >. 

— In the second, the monomer n is doubly excited; this LC is asso- 
ciated with the following local operator 

^i. = <T^l,i^"'T^na (22) 

with energy given by e^l = 4td. This LC introduces electronic cor- 
relation inside a particular double bond. We have named this LC 
D-LC. 

Notes that the LC where a monoexcitation take place in a double 
bond, let's call it S-LC, lower in energy than the D-LC, is in the (-) 
electron-hole symmetry sector; therefore, it will be a local constituent 
for excited states (Pleutin and Fave, 1998; Pleutin, 1999). 

These two LC are the analogue of the singlet along the rung for the 
two-leg ladder. 

3.2. Nearest-Neighbour-Intermonomer-Fluctuations 
(NNIF) 

With the two MC, the electrons are coupled by pairs on the double 
bonds and the resulting picture is the one of a molecular crystal in its 
ground state (Pleutin and Fave, 1998; Simpson, 1955; Mukhopadhyay, 
1995); the analogue for the two-leg ladder is given by the state | > 
represented in figure (3) with only singlets along the rungs. This crude 
picture is obviously not adapted for conjugated polymers where tt elec- 
trons are delocalized on the whole chain. As in the case of the two-leg 
ladder, fluctuations must be included in our treatment, we talk in this 
case of intermonomer fluctuations. This will be done in a somehow 
'minimal' way by considering the two following LC, extending over two 
double bonds only (Pleutin and Fave, 1998; Pleutin, 1999). 
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— In the first, one electron is transfered from one monomer to the 
nearest-neighbour one; this operation is associated with the local 
operator written as 

Cti = (23) 

with energy Cct = 2td + jV. This LC introduces some nearest- 
neighbour intermonomer charge fluctuation reproducing in a min- 
imal way the conjugaison phenomenom. We named this LC, Cti- 
LC. 

— In the second, some intermonomer spin fluctuation is introduced 
by combining two n.n. triplets into a singlet; the corresponding 

operator is 

(24) 

with energy given by eu = 4:td — {U — V); this LC is named TT-LC. 

There is more LC extended over two n.n. double-bonds within the 

(+) electron-hole symmetry class - and, of course, there is much more 
LC more extended - but we believe the two selected ones are the more 
important. For instance, we saw that the S-LC is in the (-) class of 
symmetry, however, with the simple product form of the electron-hole 
symmetry operator, a LC with two S'-LC - or more generally with an 
even number of S-LC - is in the (+) class of symmetry and should 
be considered here. However, this kind of LC are sufficiently high in 
energy to be reasonably neglected. 

These two LC are the analogue of the resonon configuration for the 
two-leg ladder. 

3.3. Recurrence relation for the ground state wave 

function 

With our choice of four LC, all possible electronic configurations are 
then build up. They are characterized by the number of D, Cti and 
TT-LCs, rid, net and nu respectively, and by the positions of these 
different local configurations. The positions of the D, Cti and TT-LCs 
are labelled by the coordinates z{k) {k = l,..,nd), y{j) {j = l,..,nct) 
and x{i) {i = l,..,ntt) respectively; implicitly, the necessary non- 
overlapping condition between LC is supposed to be fulfilled all along 
the paper. The electronic configurations are then expressed as 
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I x{l), x{ntt),y{l), y{nct),z{l), z{nd) >= 

(25) 

nr:\n-fini:iirr\.)C4)^V,)io> 

We look for a ground state of the following form 



rid 

I '^'^ >- ^d^iF ^ct Utt 11 ^ -'xW^V:/) I ^ 

{nd,nct,ntt} k=l 

(26) 

where Uf, Ud, ujf, Uct and uu are parameters to be determined and 
where the summation runs over all the possible configurations. 

This state is a kind of generalization of the Dimer-RVB state of the 
two-leg ladder. The singlet on the rung is replaced by a two component 
MC, which contain F-LC and D-hC, and the resonon is replaced by a 
two component NNIF, which contains Cti-LC and TT-LC. By analogy 
with the previous case, we can write a recurrent relation followed by 
the wave function (26); it is given by 

\N >= Uf\N - 1 > (g)|Fjv > +Ud\N - 1 > i^\Dn > + 

(27) 

uiF{uct\N - 2 > ^\CtN-i,N > +utt\N - 2 > 0\TTn-i,n >) 

where |F/v > and \D]\f > are the MC created by and 1)]^ respec- 
tively and \CtN-i,N > and |Trjv-i,Ar > are the NNIF created by the 
Cijv_i and TTj^_^ respectively. There are two constraints due to nor- 
malization to be fulfilled by the coefficients, explicitely u'j + u'^ = 1 and 
""ct + ""it ~ ^- Among the five parameters, three remain free and have 
to be determined variationally by minimalization of the total energy 
per unit cell. 

3.4. Recurrence relations for the energy 

The energy per unit cell for a polymer is given again by 

eoo = lim ejv = lim — — ^ (28) 

where 

En=<N\H\N> and Zn =< N\N > (29) 

These two quantities, the expectation value and the norm of the wave 
function, are solutions of a set of coupled recurrent relations, as for the 
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Tabic III. Energies per site from DMRG and 
RVA calculations, e^^'^ and respectively, 
for U = 3tandV = 1.2t 



A 






0.1 


0.330670 


0.396849 


0.3 


0.284416 


0.327386 


0.5 


0.224344 


0.254182 


1.5 


-0.175762 


-0.164041 



two-leg ladder case. In this case, they are simpler because the selected 
LC, MC and NNIF, are orthogonal; we don't have then to consider 
quantities as Djv and Ijv and, finally, we get only two coupled recurrent 
relations instead of four in the previous section 

Zn = + ujpZi\f_2 (30) 

En = En-1 + u]pEN-2 + f-rri{uf)ZN-l + W2{u f ,UiF ,Uct)Z ^ -2 (31) 

where €m{uf) is the correlated ground state energy of a double bond 
given by 

em{uf) = -Aultd + UfUd{U - V) (32) 

and W{uf,uiF,Uct) expressed the energy of the NNIF-LC and their 
coupling with the MC 

W2{uf,UlF, Uct) = ufpul^ect + ujpU^ett + ujpUttUctVStg + 2UfUiFUctts 

(33) 

With the two recurrent relations, wc associated the following natural 
initial conditions Zq = Zi = 1 and Eq = 0, Ei = emiup)- 

The set of equations (30) and (31) together with the definition 
(32), (33) and the initial conditions aforementioned, can be solved 
analyticaly and are also really natural to implement on a computer. 

In the following, comparisons are made with Density Matrix Renor- 
malization Group results (White, 1993; Jeckelmann, 1998; Peschel, 
1999). A finite system DMRG algorithm is used to calculate the ground 
state energy for several chain lengths up to four hundred sites and the 
energy per site is then extrapolated to the infinite lattice limit. The 
numerical error in the energy per site of the infinite lattice is estimated 
to be of the order of 10"*^ or smaller. We compare the RVA results with 
these extrapolated values. 

A reasonable choice of parameters for the polyacctylcnc is given 
hy U = 3t and V = l.2t (Boman and Bursih, 1998; Shuai, 1998). In 
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table III, we show the energy per site obtained from DMRG and RVA 
calculations for different values of the dimerization parameter A. For 
A = 2, the chain is fully dimcrized; the results given by the ansatz 
(26) is then very accurate per construction. Next, the accuracy of the 
RVA results decreases when the dimerization decreases; this appear 
clearly in the table III. For values relevant for conjugated polymers, 
A larger than 0.25, the energy approaches 90% of the DMRG energy; 
we think these results satisfactory according to the relative simplicity 
of our proposed trial wave function (26) : only four local configurations 
arc retained! This wave function is, in a sense, an extrapolation of 
the wave function analysed in small cluster calculations (Chandross, 
1999), to an infinite lattice. Of course some improvements are suitable 
and may include more extended LC; by the way, works are currently 
in progress in that direction by using again the Rcciurcnt Variational 
Approach. Nevertheless, we believe the wave function (26) sufficient to 
get a first explanation of some physical phenomena as the very strong 
linear absorption due to excitonic states observed in polydiacetylenes 
(Pleutin and Fave, 1998). 

Last, it is well known that the ground state of the EPH model shows 
a Bond Order Wave if U > IV and a Charge Density Wave in the 
other case (Mazumdar and Campbell, 1985). In figure (10), we show 
the relative error of the RVA results compared to the DMRG ones for 
several choices of Coulombic parameters in function of the dimerization 
parameter A, in the two different regimes. Once again, the results 
show clearly that our ansatz is better when A increases. Moreover, 
its behaviour seems different in the two sides of the transition, this 
is natural since our ansatz is not appropriated to describe properly a 
Charge Density Wave. 



4. Conclusions and Perspectives 

To conclude, we have presented in this note the Recurrent Variational 

Approach and, more specifically, its first application to the electronic 
structure of conjugated polymers. The ansatz built here (26), shows 
some nice similarities with the trial state (3) made for the two-leg 
ladder. 

By comparisons with Density Matrix Renormalization Group cal- 
culations, we think this trial wave function encouraging and a good 
starting point for further improvements. Moreover, we believe such 
kind of simple wave functions useful to get some analytical insight into 
physical phenomena (Pleutin and Fave, 1998; Pleutin, 1999). 
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Figure 10. Relative errors of the RVA results compared to the DMRG ones. The 
unfilled symbols are for calculations where U > 2V; these choices give a BOW 
ground state. On the contrary, the filled triangles correspond to calculations where 
U < 2V; then a CDW ground state is stabilized. 



Two natural extensions of this work are currently in progress: first, 
it is possible, in some extents, to improve the simple ansatz (26) by 
including more local configurations; second, it is also possible to apply 
the method to other conjugated polymers of current interests as the 
poly-paraphenylene and the poly-paraphenylenevinylene. 
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